Maxwell's Demon at work: 
Two types of Bose condensate fluctuations 
in power-law traps 



Siegfried Grossmann and Martin Holthaus 

Fachbereich Physik der Philipps-Universitdt, Renthof 6, D-35032 Marburg, 

Germany 

holthaus@stat.physik.uni-marburg.de 

Abstract: After discussing the key idea underlying the Maxwell's 
Demon ensemble, we employ this idea for calculating fluctuations of 
ideal Bose gas condensates in traps with power-law single-particle en- 
ergy spectra. Two essentially different cases have to be distinguished. If 
the heat capacity remains continuous at the condensation point in the 
large-iV-limit, the fluctuations of the number of condensate particles 
vanish linearly with temperature, independent of the trap characteris- 
tics. If the heat capacity becomes discontinuous, the fluctuations vanish 
algebraically with temperature, with an exponent determined by the 
trap. Our results are based on an integral representation that yields 
the solution to both the canonical and the microcanonical fluctuation 
problem in a singularly transparent manner. 
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According to grand canonical statistics, the root-mean-square fluctuations SN U of the 
occupation numbers N v of an ideal Bose gases's v-th single-particle state are given 

by 11 

{8N v f = N U {N U + 1) . (1) 

This expression follows without any approximation from the grand canonical approach, 
but it faces a severe problem when applied to the fluctuations 6 No of the ground state 
occupation number N of isolated Bose gases: if the temperature approaches zero, all N 
particles of an isolated Bose gas occupy the ground state, so that the actual fluctuations 
vanish, whereas Eq. (|l|) predicts fluctuations 8 No of the order N. This seems to be one 
of the most important examples where the different statistical ensembles can not be re- 
garded as equivalent. When computing low-temperature fluctuations of the ground state 
occupation number for isolated Bose gases, one therefore has to give up the convenient 
grand canonical point of view, and to resort to a microcanonical treatment. 

Although this problem had been well recognized and discussed some time 
ago H^l, tools for computing the microcanonical fluctuations 5 No have been devel- 
oped only recently Q , spurred by the progress in preparing Bose-Einstein condensates 
of alkali atoms in magnetic traps A particularly instructive model system for 

illustrating the microcanonical approach to fluctuations 8 No is provided by N ideal 
Bosons trapped in a one-dimensional harmonic potential. Since quasi one-dimensional 
harmonic trapping potentials can be realized as limiting cases of strongly anisotropic 
three-dimensional traps P]l0|, this system is not merely of academic interest. The value 
of the model lies in the fact that it allows one to map the problem of evaluating the 
microcanonical statistics to problems also arising in analytic number theory because 
the number of microstates accessible at some excitation energy E equals the number of 
partitions of the integer n = E/(hu>) into no more than N summands, with u being the 
oscillator frequency. Using the appropriate asymptotic formulae from partition theory 
one finds that the microcanonical fluctuations SNo for this model system vanish linearly 
with temperature T 



(2) 



where ks is the Boltzmann constant, and T^f 1 denotes the temperature below which the 
ground state occupation becomes significant. As illustrated in Fig. 1, which compares 
the relative microcanonical fluctuations SNq/N to the corresponding grand canonical 
fluctuations and to the approximation (^]), for N = 10 6 particles, this approximation 
is quite good indeed. The very same result (||) has also been obtained by Wilkens (l^] 
within a canonical approach, that is, for a trap in contact with a heat bath. 
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1. Full red line: microcanonical fluctuations SNo/N for a system of TV = 10 6 ideal Bose particles 

(i) 



trapped by a one-dimensional harmonic potential [0|. The temperature To = Tq ±; denotes the charac- 
teristic temperature below which the ground state occupation becomes significant, see Eq. (JiJ). Black 
short-dashed line: grand canonical fluctuations for the same system. Blue dashed line: low-temperature 
approximation provided by Eq. (^). 

How can one generalize this findingto other trap types? A rather interesting 
suggestion has been made by Navez et al. H|. Denoting, for an ideal ./V-particle Bose 
gas in some arbitrary trap, the number of microstates with exactly N ex excited particles 
as <&(N ex \E), so that the total number of microstates accessible at the given energy E 
reads Q(E\N) — _ <fr(N ex \E), these authors consider the generating function 



T(z,E) 



E 

AT 8x =0 



(3) 



This function involves &(N CX \E) even for N cx > N, which appears to be unphysical: 
after all, the excitation energy E can not be distributed over more than the N particles. 
However, provided the microcanonical distributions for finding iV ox out of N particles 
in an excited trap state, 



p m {N m \E) 



$(N CX \E) 

n{E\N) 



N cx = 1,2,. 



(4) 



are strongly peaked around some value N cx <C N, which should be the case 
for temperatures well below the onset of Bose-Einstein condensation, we will have 
$(N ex \E) /fl(E\N) « for A^ cx > N. In that case the generating function (||) would be 
quite useful, since one could obtain the microcanonical expectation value (N ) for the 
ground state occupation number, and its fluctuation, from 



N-(N ) = z^\nT{z,E) 
oz 



and (SN ) 2 = ( Z JL) lnT(z,E) 
z=i V dz 



z=l 



respectively 14 . The proviso can be formulated in more intuitive terms: the required 
well-peakedness of the distributions (^|) means that those microstates where the energy 
E is actually spread out over all N particles carry only negligible statistical weight, 
so that the overwhelming majority of all microstates leaves a fraction of the particles 
in the ground state, forming the Bose condensate. Then the restriction on the number 
of microstates caused by the fact that there is only a finite number N of particles 
becomes meaningless, so that, loosely speaking, "the system has no chance to know 
how many particles the condensate consists of". But if this is the case, i.e., if the 
system's properties become insensitive to the actual number of particles contained in 
the condensate, then one can act as if the condensate particles constituted an infinite 
reservoir. Thus, the generating function (|^) may be regarded as the partition function 
of a rather unusual ensemble, consisting of the excited-states subsystems of Bose gases 
that exchange particles with the ground state "reservoirs" without exchanging energy. 
Since such an exchange process, if performed by hand, requires a genius who is able 
to separate the hot, excited particles from the cold ones in the ground state, this new 
ensemble has been called the "Maxwell's Demon ensemble" [H. 

But can we rely on Maxwell's Demon, that is, does the proviso hold? This 
question needs to be answered first. A strong argument in favour of the Maxwell's Demon 
ensemble has already been provided by the approximation (0) to the low-temperature 
fluctuations for a Bose gas in a one-dimensional oscillator potential: these fluctuations 
are independent of the total particle number N, as they should be if the system really 
has no knowledge of the number of condensate particles, and thus of N. The awe- 
inspiring agreement with the actual microcanonical fluctuations depicted in Fig. 1 leaves 
no doubt that this approximation is reliable. To further substantiate the new ensemble, 
we also consider the microcanonical fluctuations 5 No for an ideal Bose gas trapped by a 
three-dimensional isotropic harmonic oscillator potential |15| ]. The numbers Q(E\N) of 
microstates for some given excitation energy E = nfwj can then be obtained from the 
canonical A^-particle partition function 

oo 

Zn{P) = 53 e ~ n>3n " n(nKu\N) , (5) 

n=0 

which, in turn, can be calculated numerically with the help of the recursion rela- 
tion 



1 N 

znW) = ^j2 Zi ^ ZN - k ^ ■ ( 6 ) 



fc=l 



As usual, (3 = 1/(&bT) denotes the inverse temperature. By means of numerical saddle- 
point inversions of Eq. (^|), we compute the desired numbers Sl(E\N cx ) for N cx ranging 
from 1 to N Jig ], and get the differences 

$(iV ex |£) = Sl(E\N cx ) - fl{E\N cx - 1) (7) 

that determine the microcanonical distributions (^). Some of these distributions are 
displayed in Fig. 2, for N = 1000 and several "low" temperatures. What we find is 
exactly what is needed for Maxwell's Demon: the distributions are well peaked for 
temperatures below the onset of condensation, and remarkably close to Gaussians jl2| . It 
is then no surprise that the corresponding microcanonical low-temperature fluctuations 
SNq, obtained from the widths of these distributions, are — as long as a condensate 



exists! — once again independent of N, as exemplified in Fig. 3 for N = 200, 500, and 
1000. As discussed above, it is precisely this iV-independence, expressed mathematically 
by the appearance of the upper summation bound "oo" rather than "iV" in Eq. (|3|), that 
lies at the bottom of the Maxwell's Demon ensemble. But whereas this iV-independence 
is, by construction, put into this ensemble, it has come out here as the result of a truly 
microcanonical calculation [|l2|Jl5f| that works with the actual N, not with oo. 




N ex /N 

Fig. 2. Microcanonical probability distributions p C x(Af ex |n) for finding N ex out of JV = 1000 ideal Bose 
particles, trapped by a three-dimensional isotropic harmonic potential, excited when the total excitation 
energy E is n ■ hw, with ui denoting the oscillator frequency. The number n determines the temperature 
T. The normalized temperatures T/Tq corresponding to the blue, Gaussian-like distributions range 
from 0.3 to 0.9 (left to right, in steps of 0.1); T = (hu>/k B )(N/C(3)) 1 ^ 3 . Due to finite- JV-effects, the 
condensation temperature is lowered from To to about 0.93 To . The temperature corresponding to the 
rightmost, red distribution is T = 0.95 To, lying slightly above the condensation point. 
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Fig. 3. Microcanonical fluctuations SNo for iV = 200, 500, and 1000 ideal Bose particles trapped by a 
three-dimensional, isotropic harmonic potential. The fluctuations are maximal close to the respective 
condensation points. These maximal fluctuations scale approximately as \f~N, cf. Eqs. (^) and (^fj). 
Note that the low-temperature fluctuations for all three systems agree perfectly, thus demonstrating 
the ^-independence of (STVq below the condensation point. 



Having thus gained confidence in the abilities of Maxwell's Demon, we now set 
it to work in order to compute condensate fluctuations 5N . To this end, we consider 
ideal Bose gases in d-dimensional traps with arbitrary single-particle energies £„; we 
stipulate Eq — 0. Denoting the grand canonical partition function by 3(2, 0), we base 
our analysis on its "excited" part S ex (^, (3) = (1 — z)E(z, /?). Since, by virtue of Eq. (|7|), 

oo 

S cx (z,/3) = J2 z N »*J2$(N ex \E)e-P E , (8) 

7V QX =0 E 

this function has the decisive property 



. * \ k 



z-q- z ) 3 ex (z,/3) 



E e-(> E = MfcOS) , (9) 



E \JV cx =0 



i.e., it yields directly the non- normalized canonical moments Mk{0), and generates the 
microcanonical moments 

oo 

» k (E) = J2 N^(N^\E) with fc = 0,l,2,.... (10) 

N ox =0 

We then employ the Maxwell's Demon approximation: as long as there is a condensate, 
these moments (with iV 0X ranging from to oo) approximate the true moments of the 
physical set { Q(N CX \E) } (where the number iV cx of excited particles can not exceed 
the total particle number N); in this approximation one has the identity ^o{E) — 
Q(E\N) Jig] , Now the calculations within the canonical ensemble become remarkably 
simple. The canonical expectation value (N cx ) = N — (No) of the number of excited 
particles is given by 

MAQ) , s 

<"«> = AW) ; 

the canonical condensate fluctuations (5N )l n — {5N cx ) 2 cn follow from 

{6N0) ™ = MoM " ImoMJ ' (12) 

Without any further approximation, these expressions can be rewritten as complex 
integrals: 

i />T+ZOO 

(iVex) = 7T- / dtT(t)Z((3,t)C(t) (13) 

and 

-i />T+ZOC 

(M f cn = — dtr(t)Z((3,t)C(t-l) , (14) 

2.m J T _ loo 

where T(t) and £(i) denote the Gamma function and Riemann's Zeta function, respec- 
tively. All the information about the specific trap under consideration is embodied in 
its spectral Zeta function 

oo 1 

mt) = d5) 



where the sum runs over the trap spectrum, excluding the ground state energy Eq = 0. 
The real number r in Eqs. ( |l3|) and (|l4| ) has to be chosen such that the path of inte- 
gration up the complex i-plane sees all poles to its left. 

So far, the analysis is quite general. We now specialize the further deliberations 
to ideal Bose gases in d-dimensional traps with power-law single-particle spectra 



A^c^f , Vi = 0,1,2,... , a >0 , (16) 



where the dimensionless coefficients Cj characterize the trap's anisotropy, normalized 
such that the lowest Ci is unity; the characteristic energy A measures the gap between 
the ground state and the first excited state, and the exponent a is determined by the 
potential's shape. Such systems have been studied first by de Groot et al. Jl9|] ; we have 
adopted here the notation also employed by Wilkens and Weiss (TsJ . 

If we consider 7V-asymptotically large systems and disregard finite- TV-effects, 
that is, if we focus on gases consisting of at least some 10 5 particles, say, then a good 
approximation to the density of states is provided by 

r(|) W (nt,*) 

Using this densit y, a nd assuming that the anisotropy coefficients Cj are not too different 
from each other |10[ , the usual line of reasoning shows that for d/a > 1 there is a sharp 
onset of Bose-Einstein condensation at the temperature To given by JTTj[] 



Moreover, the spectral Zeta functions can now be well approximated by 



(18) 



Z(P,t) ^ ^-(^Ay'at+l-d/a) , (19) 



so Eqs. (|l^ ) and (14) adopt the transparent forms 



A 



^ ~ f(l)2^ / ^ (/3A)-*r(t)C(* + l-rf/a)CW (20) 



and 



4 1 rT+ioo 

(^o)L w fviy ^- y d * ^ A r ( r w c(* + 1 - d/^) c(* - 1) • (2i) 

For A -c fc^T, the behavior of either integral is determined by the pole of its 
integrand farthest to the right in the complex plane. Keeping in mind that (^(z) has 
merely one single pole at z = 1, with residue 1, while the poles of T(z) are located at 
z = 0, —1, —2, . . . , the decisive pole is provided either by the system's Zeta function 
((t + 1 — d/a), or by the other Zeta function that is determined by the order of the 
cumulant one is asking for: by ((t), if one asks for the first cumulant (N cx ), or by C(* — l)i 
if one asks for the second cumulant (SNp) 2 . To see what the argument boils down to, let 
us first consider the evaluation of Eq. (|20|), where the system's pole at t = d/a competes 
with the cumulant-order pole at t = 1: 



• If d/a > 1, the low-temperature behavior of (N ex ) is governed by the pole of 
((t + 1 — d/a) at t — d/a. Hence the residue theorem yields 



(iVe: 



A ((d/a) 



k B T 



d/a 



(22) 



This canonical result, valid for T < Tq, coincides precisely with the result of the 
customary grand canonical analysis p9| . For example, in the case of the three- 
dimensional isotropic harmonic oscillator potential (i.e., for d = 3, a = 1, A = 1 
and A = Tilj) Eq. (E2) yields the familiar formula 



(N ) = N - (N m ) 



N 



1 



for T<T 



Huj ( N 



C(3) 



1/3 



• If d/a = 1, both Zeta functions in Eq. (|20|) coincide. We then encounter a double 
pole at t = 1, and find 



(iVe: 



.4 



k B T 
A 



(23) 



where 7 = 0.5772 ... is Euler's constant. This corresponds to a result obtained 
already in 1950 by Nanda |2(J with the help of the Euler-Maclaurin summation 
formula. 



If < d/a < 1, the pole of £(£) at t = 1 takes over: 



A 



r(f) 



C(2 - d/a) 



(24) 



so that for sufficiently low temperatures (iV ex ) now depends linearly on T, regard- 
less of the value of d/a that characterizes the trap. 

A mere glance at Eq. ( pl| ) then suffices to reveal that the very same scenario — a first 
pole at t = d/a that endows the temperature dependence with a trap-specific exponent 
as long as it lies to the right of a second one, which yields universal behavior when it 
becomes dominant — also governs the canonical condensate fluctuations, with the only 
difference that the second pole now is located at t = 2: 

• If d/a > 2, the pole of £(£ + 1 — d/a) at t — d/a wins, giving 



(SNo)i 



AQ(d/a-l) 



d/a 



(25) 



If d/a = 2, we find at t = 2 the already familiar double pole, resulting in 

2 



(SN t 



0) 



A 



A 



In 



k B T 
A 



+ 7 + 1 



(26) 



• If < d/a < 2, the pole of £(i — 1) at t = 2 lies to the right of its rival, yielding 

In particular, for the one-dimensional harmonic oscillator we have d = 1, a = 1, 
A = 1 and A = hui, so that we recover our previous microcanonical result (|^) 
within the canonical ensemble, recalling that £(2) — ir 2 /6. 



The above canonical fluctuations, derived from the integral representation (|i"4|), reduce 
to the expression obtained by Politzer |2l[] in the case of the three-dimensional isotropic 
trap, and match the results obtained by Wilkens and Weiss |l8| . 

The calculation of the corresponding microcanonical quantities now requires 
saddle-point inversions of Eq. (Q) in order to obtain the microcanonical moments /i^ (E) 
from the canonical moments Mk(f3). Performing these inversions, and reexpressing en- 
ergy in terms of temperature, we find that the integral (jl^) — and, hence, the re- 



sults (|22|), (23), and ( |24| ) for the number of excited particles — remains valid within 
the microcanonical ensemble. The fluctuations require more care: Whereas canonical 
and microcanonical fluctuations coincide in the large- iV-limit for d/a < 2, the micro- 
canonical mean-square fluctuations (SNo) mc are distinctly lower than their canonical 
counterparts for d/a > 2: 



Ad C 2 (f) (k B T 



(5N y cn - (5N y mc « — ^ffi [-^~ ) - ? ^'"1 T T { , 



d/a 



(28) 

Thus, the exponent of T is the same for both (5N ) 2 n and {5N ) 2 mcl but the prefactors 



can differ substantially. This Eq. ([28|) contains as a special case the result obtained for 
the three-dimensional isotropic trap by Navez et al. p3|. 

Before summarizing these findings, it is useful to also consider the heat capaci- 
ties for trapped ideal Bose gases with the single-particle spectra (|l6|): for d/a > 1, and 
temperatures below the condensation temperature T , the heat capacity per particle is 
given by 



Nk B a\a J C(|) V^o 
above T by 

-% = + g ^ + ; ( f-4 9d/Az \- (30) 

Nk B a\a J g d / a (z) a 1 g d / a -i{z) 

Since the fugacity z approaches unity from below when T approaches To from above, 
so that the Bose function g a {z) approaches C( a )> we see that the heat capacity remains 
continuous at Tq for < d/a < 2, but exhibits a jump of size 



Nk B 



c(-) 

(31) 



To 



o* C (2-l 



for d/a > 2. 

We thus arrive at the following picture: For any dimension d and trap expo- 
nent a > 0, the fluctuation of the number of condensate particles is independent of 
the total particle number N. For isolated traps, this insensitivity of the system with 
respect to N reflects the well-peakedness of the microcanonical distributions (Q), see 
Fig. 2: if there is a condensate, the behavior of the ideal Bose gas does not depend on 
how many particles the condensate consists of. If d/a < 2, so that the heat capacity 
remains continuous in the large- AT-limit, canonical and microcanonical fluctuations 5 No 
vanish linearly with temperature, see Eq. (|27j). If d/a = 2, there appears a logarithmic 
correction to the linear T-dependence, as quantified by Eq. (p6|). But if d/a > 2, so that 
the heat capacity becomes discontinuous, then the fluctuations <5iVo vanish proportion- 
ally to T d / 2CT , so that now the properties of the trap determine the way the fluctuations 
depend on temperature. In addition, in this case the microcanonical fluctuations are 
markedly lower than the fluctuations in a trap that exchanges energy with a heat bath. 



Intuitively, one might have expected some sort of square root law for the fluctu- 
ations. Because of the iV-independence of the condensate fluctuations, there is, of course, 
no "v^/V-dependence" of SNq. The square root is hidden elsewhere: Since SNq = 5N ex 
for ensembles with fixed particle number N, we find 

• 6N CX cx (7V C x) 1/2 for 2 < d/a; 

• 5N CX cx (N cx y /d for 1 < d/a < 2; 

• 5N CX oc (N cx ) for < d/a < I, 

with proportionality constants that are independent of temperature, both canonically 
and microcanonically. The first of these relations is just what one might have guessed, 
but the crossover from normal fluctuations for d/a > 2 to much stronger fluctuations 
for < d/a < 1 appears noteworthy. 

It remains to be seen how much of this ideal structure survives in the case 
of weakly interacting Bose gases. It should also be recognized that Maxwell's Demon, 
though it has provided the microcanonical low-temperature fluctuations, can not solve 
all problems of the ideal gas. When considering <i-dimensional isotropic harmonic traps, 
the Maxwell's Demon approximation (i.e., the replacement of the true upper summation 
bound "N" in Eq. (|l0|) by "oo") is exact below the "restriction temperature" (i.e., that 
temperature where the number n — E/(huj) of energy quanta equals the number N of 
particles |l5|]), but the description of the Bose-Einstein transition itself is beyond the 
capabilities of Maxwell's Demon. Namely, that description requires the computation 
of the numbers £l(nh~Lu\N) of microstates also under conditions where the restriction 
due to the finite N becomes decisive. Incidentially, one meets the task of computing 
such restricted partitions of integers also in other problems of statistical mechanics, for 
example in the theory of the so-called compact lattice animals, or of the infinite-state 
Potts model 

Nonetheless, the results obtained with the help of the Maxwell's Demon approx- 
imation have some interesting number-theoretical implications. Going once more back to 
Eq. (0) for the one-dimensional oscillator, and inserting the energy-temperature relation 
n = E/(huo) sa £(2)(fcsT/?kj) 2 , we find the truly remarkable formula 

5N a « ^fh • (32) 

This has a twofold interpretation. The physicist, puzzeled by the loss of the square root 
fluctuation law at the level of (N ex ), finds a substitute: 

• For ideal Bose particles trapped at low temperatures by a one- dimensional har- 
monic potential, the root-mean-square fluctuation of the number of ground state 
particles is given by the square root of the number of energy quanta. 

The mathematician, who approaches Eq. (|3^) from the viewpoint of partition theory, 
sees the solution to another problem: 

• If one considers all unrestricted partitions of the integer n into positive, integer 
summands, and asks for the root-mean-square fluctuation of the number of sum- 
mands, then the answer is (asymptotically) just y/n 

— certainly one of the most amazing examples for the occurrence of square root fluc- 
tuations! The ease with which the solution to a seemingly difficult number-theoretical 
question has been obtained here is even aesthetically appealing. It is pleasing to conclude 
that ongoing developments in statistical mechanics, themselves being motivated by re- 
cent experimental achievements ^,^,||, have a high potential for further fertilization 
across subfield boundaries. 



